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Abstract 

In this second paper on quantum fluctuations near the classical instanton configurations, see |l|, 
we focus on another well studied quantum-mechanical problem, the one-dimensional Sine-Gordon 
potential (the Mathieu potential). Using only the tools from quantum field theory, the Feynman 
diagrams in the instanton background, we calculate the tunneling amplitude (the instanton den¬ 
sity) to the three-loop order. The result confirms (to seven significant figures) the one recently 
recalculated by G. V. Dunne and M. Unsal, Phys. Rev. D 89, 105009 (2014) from the resurgence 
perspective. As in the double well potential case, we found that the largest contribution is given 
by the diagrams originating from the Jacobian. We again observe that in the three-loop case 
individual Feynman diagrams contain irrational contributions, while their sum does not. 
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Introduction 


Since it is our second paper of the series, following the one on the double-well potential 
|l| , it does not need an extensive introduction. Topological solitons, instantons in particular, 
are widely used in the context of quantum held theories and condense matter physics. Their 
relation to standard perturbative series is an old issue, which continues to produce interesting 
results, so far mostly in quantum mechanical context. 

The Sine-Gordon (SG) held theory has been extensively studied in classical context, 
with an enormous literature dedicated to it, see e.g. and references therein. Goleman 
has extended the results to the quantized theory by relating the SG held theory to the 
zero-charge sector of the massive Thirring model. Also in [d] the explicit calculations for 
the tunneling amplitude using the so-called nonvacuum instantons at hnite energy were 
presented. 

The quantum mechanical SG potential (the Mathieu potential) is the basic element of 
condense matter theory. Tunneling from one minimum to the next, in the path integral 
formulation, is described by Euclidean classical paths - the instantons. The issues we dis¬ 
cuss in this paper deal with quantum huctuations around these paths. We would like to 
demonstrate by an explicit calculation how our tools work in this - well controlled and 
studied setting -before applying them to more complicated/realistic settings in quantum 
held theory. Therefore we do not use anything stemming from the Schrodinger equation in 
this work, in particularly do not use series resulting from recurrence relations or resurgence 
relations (in general, conjectured) by several authors. 

One reason to study SG is to explore further the existing deep connections between 
the quantum mechanical instantons - via Schrodinger equation - with wider mathematical 
issues, of approximate solutions to diherential equations, dehned in terms of certain gen¬ 
eralized series. A particular form of an exact quantization condition was conjectured by 


B. 


J. Zinn-Justin |5|1, which links series around the instanton, instanton/anti-instanton sectors 
with the usual perturbative series in the perturbative vacuum. It remains unknown whether 
it can or cannot be generalized to the held theory cases we are mainly interested in. Re¬ 
cently, for the quartic double well and Sine-Gordon potentials Dunne and Unsal (see jol 
and also references therein) have presented more arguments for this connection and made it 
more precise, which they call resurgent relation between perturbative and instanton sectors. 
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Another reason for which we decided to do this work is a certain set of observations abont 
Feynman diagrams on top of the instanton for the double well potential with degenerate 
minima we observed in our first paper [l|. We wanted to see how general they are, using 
another example, now with infinitely many degenerate minima. The SG potential also has 
new vertices and thus many new diagrams. As we will show below, indeed all these trends 
repeat themselves in this second setting as well. 

Few comments on the history of present approach. Omitting well known classic papers 
on instanton calculus we mention a pioneering paper j^, where the two-loop correction to 
the tunneling amplitude for the SG was calculated. In particular, the formalism for treating 
the zero-mode singularities was described in detail. 


Three-loop correction to the instanton density 


Let us consider the quantum-mechanical problem of a particle of mass m = 1 in the Sine- 
Gordon potential 

V = ^[1 -cos(^x)] . (1) 

The well-known instanton solution Xinst{t) = ^ arctan(e*) describes the tunneling between 
adjacent minima by the Euclidean classical path with the action Sq = S[Xinst{t)] = p- Our 
notation for the coupling is related to those used in Q by = |. The classical action Sq 
of the instanton solution is therefore large and ^ is used in the expansion. 

The SG potential has an infinite number of degenerate minima, and perturbative levels 
in them form a continuous band, with states within the band labeled by Bloch angular 
parameter 6. The energy of the lowest band is 


^(lowest band) _ , 


6E 

T 


( 2 ) 


where Eq is the naive ground state energy, without tunneling, written as the following 
expansion 

- OO ^ 

Eo = = P) 


2 ^ Si 

n=0 U 


1 .1 r 7-T 7 -i(lowest band) 7 -i(lowest band) , , i . i j. i 111 

while 0 E = Eg^^ — Eg^Q generates another series, related to the so called 
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instanton density 


OO 


6E = AEJ2 


On 
-n *^0 


{Bo = 1 ) , 


( 4 ) 


n=0 


here AE = 2-^^ e is the well-known one-loop semi classical result [^. Coefficients 
in the series (jS]) can be calculated using the ordinary perturbation theory (see 3) while 
many coefficients in the expansion (0]) were found by J. Zinn-Justin and collaborators, 
1981-2005 (see 3 references therein), obtained via the so called exact Bohr-Sommerfeld 
quantization condition. It was remarkably refined by Dunne-Unsal 3 ^^e formalism of 

the so-called resurgent trans-series. 


Alternatively, using the Feynman diagrams technique Lowe and ^one 3 calculated the 


two-loop correction Bi = —7/8 which was later on reproduced in [5[ in the so-called exact 
Bohr-Sommerfeld quantization technique. Higher order coefficients in (j4]) can also be 
computed in this way. Since we calculate the energy difference, all Feynman diagrams in the 
instanton background (with the instanton-based vertices and the Green’s function) need to 
)e accompanied by subtraction of the same diagrams for the trivial x = 0 saddle point (see 
3 for details). For r 3> 1 it permits to evaluate the ratio 


(7r|e 

(0|e-^lO),=o 

where the matrix elements > ( 0 |e“^'’'| 0 ) 3 ;=o are calculated using the 

instanton-based Green’s function and the Green function of the harmonic oscillator, re¬ 
spectively. 

The instanton-based Green’s function G{x, y) form to be used 


G{x,y) = 


G^{x,y) 


1 4 X y x"^ y^ x"^ y 


2(1 + x2)(l + 2/2) _ 


+ (1 — Axy + x^y^ + + 2(1 — xy)\x — y\) log(2 G^(x, y)) 


( 5 ) 


is expressed in variables x = tanh(Y), y = tanh(y) , in which the familiar Green function 
ig-|ti-t 2 | Qf harmonic oscillator is 


G%x,y) 


1 — \x — y \ — X y 
2(1 + \x — y \ — X y) ’ 


( 6 ) 


In its derivation there were two steps. One was to find a function which satisfies the Green 
function equation, used via two linearly-independent solutions and standard Wronskian 
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method. The second step is related to a zero mode: one can add a term 0 o(^i) 0 o(^ 2 ) with 
any coefficient and still satisfy the equation. The coefficient is then fixed from orthogonality 
to the zero mode. 

The two-loop coefficient in (|1]) is 


Bi — a -\- h\ -\- 1)2 -\- c ^ 
53 , 3 ,7 

CL —-, 0i — — , bo — — ^ 

60’ 40 ’ 20 ’ 


c = 


5 

12 


( 7 ) 


reflecting the contribution of four Feynman diagrams, see Fig. [T] 


The three-loop correction B 2 (S]) we are interested in is given by the sum of twenty-two 
3-loop Feynman diagrams, which we group as follows (see Figs. [2]-[3]) 

B2 = Oil bii bi 2 -|- 621 + ^22 + ^23 + ^24 

( 8 ) 

Cl-\-C2-\-c^-\-€4 ^-\-c^-\-Cq-\-Cj-\- B2100P > 
complementing by a contribution from two-loop Feynman diagrams, see Fig. [T|, 

1 , , r \2 / 7 7 \ 341 

B 2100 P = —(0 + O 1 -I-O 2 ) -l-(o-l-Oi-l- 02 ) c = ’ 

(see (Ej)). 

The rules of constructing the integrals for each diagram should be clear from an example: 
the explicit expression for the Feynman integral 623 in Fig. |2l which is 


623 = 32768 / dx j dy 


dz 


dw J{x,y,z,w) [xyzw 


while for C 4 in Fig. [3] it takes the form 


>1 




C 4 = 256 / dx I dy 


dz- 


xy{l — 62 ;^ - 1 - z^) 


j-I j-I 

here we introduced notations G 


_i {l + x^f{l + y^f{l + z^f{l-z‘^) 


G (^2 

xy ^yz ^ ZZ 5 


( 9 ) 


( 10 ) 


xy 


= G{x,y),Gly = G°{x,y) and J = 

^■ysi Notice that the c's diagrams come from the Jacobian of the 
{i+x^Y {i+y^Y Y+^^Y (i+«i^) ^ 

zero mode and have no analogs in the perturbative vacuum problem. 

For calculation of the symmetry factors for a given Feynman’s diagram we use the Wick’s 
theorem and contractions, see e.g 8 |. It can be illustrated by the next two examples. For 
diagram j the three bubbles can be rearranged in 3! ways and each propagator, which starts 
and ends at the same vertex forming loop, contributes with a factor of two giving a symmetry 
factor of 2^ X 3! = 48 . For the diagram C 3 with no bubbles the three propagators which 
connect the same pair of vertices can be rearranged in 3! = 6 ways only. 
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Results 


The obtained results are summarized in Table HI All diagrams are of the form of one¬ 
dimensional, two-dimensional, three-dimensional and four-dimensional integrals. The five 
diagrams fen, d, k, I, C 7 , in particular (see Fig. |2]) 


bn = 


16 


dx / dy 


'-1 

rl 


d =16 


dx / dy 


-1 (1-a:2)(l-1/2) 

1 


(l_ 6 x 2 + a;4)(l-6|/2 + //4) 4 


(H-x 2 ) 2(1 + ^ 2)2 


■xyj 


>-1 


-1 ( 1 - 2 ; 2)(1 - 1 / 2 ) 

(1 — 6x^ -1- x'^)(l — Qy^ + y'^) 

(l + a;2)2(l +1/2)2 


0 ^ 2^0 


xy'^yy '^xxy^xy) ^ 


yy 


( 11 ) 


correspond to two-dimensional integrals and together with diagram j, a one-dimensional 
integral, are the only ones which we are able to calculate analytically 


bn = 
d = 


180100 1 , 

+- 178 C(2) - 204 C(3) + 27C(4) = b[f + 6*/™* 


756000 900 

73931 


-'ll 


47250 900 

184 

J = — , k = 


( 12 ) 


379 


I = 


244 


C 7 = 


16 


315 ’ 630 ’ 945 ’ i —1 ’ 

here ({n) denotes the Riemann zeta function of argument n (see (ill). Diagrams bn, d, 
contain a rational and an irrational contribution such that 

^irrat 


yirrat 

b\f 


-0.341 


drat 


-0.338 


It shows that for diagrams 5ii and d the rational contribution is three times larger than 
the irrational part. In the case of the DW potential the situation is opposite, the irrational 
part is dominant (see Q,. Olhe. c9ag.a.s. see Table , wese evaluated uuusencally with au 
absolute accuracy ~ 10“®. Surprisingly, almost all of them (20 diagrams out of 22 ones in 
total) are of order > 10 “^ as for B 2 itself with two of them (diagrams bn, 621 ) which are of 
order 10 “^. 

Dunne-tinsal (see jo] and references therein) reports a value of 

59 


Dunne—Unsal 


B. 


while present calculation shows that 


128 


= -0.4609375 , 


(13) 


represent 

^2 


-0.460937498 , 


(14) 
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which is in agreement, up to the precision employed in the numerical integration. 

Similarly to the two-loop correction Bi the coefficient B 2 is negative. For not-so-large 
barriers (Fq ~ 1 ), the two-loop and three-loop corrections are of the same order of magnitude. 
The dominant contribution comes from the sum of the two-vertex diagrams d, fen,/c,/, Cy 
while the four-vertex diagrams fei 2 , fe 2 i, fe 23 , e, h, Ci, C 5 , Ce provides minor contribution, the 
absolute value of their sum represents less than 0.2% of the total correction B 2 . Interestingly 
for both two and three loop cases the largest contribution comes from the ’ears’-like diagrams 
a and d, respectively, ^ 1.01 and 2.25 . 

We already noted that individual three-loop diagrams contain irrational numbers. Since 
the Dunne-Unsal’s result is a rational number, then there must be a cancelation of these 
irrational contributions in the sum ([H]). From (IT^ we note that the term gives a 

contribution of order one to the mentioned sum (| 8 ]), and therefore the coincidence in the order 
of 10“® between present result flT4|) and one of Dunne-Unsal (ir3|) is an indication that such a 
cancelation occurs. Now, we evaluate the coefficients Ai, A 2 in ([3]) using Feynman diagrams 
(see [9|). In order to do it let us consider the Sine-Gordon potential V = ^[1 — cos{gx)] 
and calculate the transition amplitude {x = 0|e“'^'^|x = 0). All involved Feynman integrals 
can be evaluated analytically. In the limit r —)■ cxo the coefficients of order Sq^ and Sq'^ in 
front of r gives us the value of Ai and A 2 , respectively. As it was mentioned above the c's 
diagrams do not exist in this case. The Feynman integral a in Fig. [Ugive us the value of 
Ai, explicitly it is equal to 

Ai = -2 . 

The diagrams fen, d and j in Fig. [2] determine A2, fen = —| d = —8 and j = y. Then 

Ao = -4 , 


which is in agreement with the results obtained in standard multiplicative perturbation 
theory (see e.g. 12|). No irrational numbers appear in the evaluation of Ai and A 2 . 


Conclusions and Discussion 

In conclusion, we have calculated the tunneling amplitude (level splitting related to the 
instanton density) up to three-loops using Feynman diagrams for quantum perturbations on 
top of the instanton. Summing all of these contributions we obtain the third coefficients B 2 
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(defined in (j4])). The resnlt - to the nnmerical accnracy we kept - is fonnd to be in good 
agreement with the resurgent relation between perturbative and instanton series suggested 
by J. Zinn-Justin and collaborators, and Dunne-Unsal (for modern reference see Q|). 

Let us remind again, that this paper is methodical in nature, and its task was to develop 
tools to calculate tunneling phenomena in multidimensional QM or QFT context, in which 
any results stemming from the Schrodinger equation are not available. We use a quantum 
mechanical example as a test of the tools we use: but the tools themselves are expected to 
work in much wider context. 


Q) 


When we started these works (see llH) we, naively, expected to see some correspondence 
between vacuum and instanton series on the level of individual Feynman diagrams. However, 
no such trend has been detected so far. Furthermore, “new” diagrams originating from the 
instanton zero mode Jacobian, surprisingly, provide the signihcant contribution ~ 50% to 
the two-loop correction Bi (one diagram out of four, see Fig.l) and ~ 114% to the three-loop 
correction B 2 (seven diagrams out of 22, see Figs.2-3), see Table I, both for 2-loop and 3-loop 


contributions, both for the double well [l| and SG problems. In the double well case 
“new” individual diagrams c and C 5 give the dominant contribution (~ 83% and ~ 127%) to 
the overall loop coefficients Bi and B 2 out of 4 and out of 18 diagrams, respectively, while 
in the SG case they give signihcant contributions ~ 48% and ~ 25% out of 4 and out of 22 
diagrams, respectively. However, the corresponding cs-like four-loop diagram in the SG case 
represents the ~ 4% of the four-loop correction B 3 only. We calculated C 2 -like three-, four-, 
hve-loop diagrams (see Fig. 3 for 3-loop case and Note Added as of Nov.l, 2015, see below): 
This single tadpole diagram gives ~ 75%, ~ 50% and ~ 30% contribution for three-, four-, 
hve-loop cases, respectively. It is quite amusing that in three-loop case the sum of C 2 and G 5 
diagrams gives ~ 100%. It implies that the sum of remaining 20 diagrams is almost zero! 

Another observation is that the hnal three-loops answer has a rational value. However, 
unlike the evaluation of the two-loop coefficient Bi where all Feynman diagrams turned out 
to be rational numbers, in our case of B 2 at least two diagrams contain irrational parts. 
What is the origin of these terms and how they cancel out among themselves are questions 
left unanswered above, since several diagrams had resisted our ehorts to get the analytic 
answer, so that we used numerical multidimensional integration methods, in particular, a 


Q 


the 


dynamical partitioning 


13| . Perhaps, this can still be improved. 


Similar calculations for scalar and eventually gauge theories would be certainly possible 






and are of obvious interest. The diagrams are the same, and the basic element remains 
explicit Green functions. (In the case of gauge theories those should be orthogonal to all 
-including gauge- zero modes.) 
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Note added in proof (July, 2015): 

Note that the normalization of the Green functions G in this paper is different from 

n 

1|. After the paper was submitted we obtained additional results. Using Dunne-Unsal 
procedure we calculated the R 3 , B 4 coefficients in the expansion (jlj) 

^3 = -897/1024 

Bi = -75005/32768 

which again turned to be rational. We evaluated contributions of the c, cs-like diagrams, 
with maximal number of integrations, to the four and five loop coefficients. Those diagrams 
still contribute significantly to the total answer, although the number of diagrams grows 
dramatically with order. Those are ~ 48%, ~ 25%, ~ 4%, ~ 0.4% of total two-, three-, 
four-, hve-loop Ri, R 2 , 7 ? 3 , R 4 coefficients, respectively. Unlike the case of the double-well 
potential [l| the absolute values of these diagrams tend to decrease by factor ~ 3 with each 
order. 

Note added (November, 2015): 

We evaluated contributions of the C 2 -hke diagrams (see Fig.3), with maximal number of four- 
point vertices V4, to the four and five loop coefficients. Those are ~ 75%, ~ 50%, ~ 30% 
of total three-, four-, hve-loop B 2 ,B^,B 4 coefficients, respectively. These anomalously large 
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contributions need to be explained. 
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Figure 1: Diagrams contributing to the two-loop correction Bi = a + bi + b 2 + c. They enter into 
the coefficient B 2 via the term B 2 ioop- For the instanton field the effective triple, quartic, quintic 
and sextic vertices are I /3 = 4 ^2 F 4 = 8 - 1) Sq"' , F 5 = -32 ^2 , 

Ve = —64 — 1) , respectively, all marked by (filled) bullets, while for the subtracted 

vacuum field diagrams we have V 3 = V 5 = 0, V 4 = 8 S'(C^ and Vg = 64 . The tadpole in diagram 

c, which comes from the zero-mode Jacobian rather than from the action, is effectively represented 
by the vertex (Jacobian source) Vtad = ^ marked (unfilled) open bullet (see text). 

The signs of contributions and symmetry factors are indicated. 
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Feynman 

diagram 

Contribution to 

B 2 

Oi 

0.076497376 

bi2 

0.012018080 

b2i 

0.021795546 

b22 

0.126625453 

b23 

0.089343699 

b2A 

0.136767200 

e 

0.095677222 

f 

0.316847853 

9 

0.308988525 

h 

0.061784864 

Cl 

-0.053695603 

C2 

-0.347282176 

C 3 

-0.071612992 

C4 

-0.183024909 

C 5 

-0.114906259 

C6 

-0.111263501 

I 2 D 

-1.70563 

Hd 

0.36381 

Iad 

0.00075 


Table I: Contribution of each diagram in Fig. [2]-[3]for the three-loop correction B 2 with symmetry 
factor included. We write B 2 = {B 2100 P + ho + ho + ho + I ad) where j = ho and hn^hD, I ad 
denote the sum of two-dimensional, three-dimensional and four-dimensional integrals, respectively. 
The term B 2100 P = 341/1152 « 0.296 (see text). 
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Figure 2: Diagrams contributing to the coefficient i? 2 - Triple, quartic, quintic and sextic vertices 
V 3 ,V 4 ,V 5 ,V 65 ah are marked by (filled) bullets. The signs of contributions and symmetry factors 
are indicated. 


13 





















Figure 3: Tadpole diagrams contributing to the coefficient B 2 . They come from the Jacobian of 
the zero mode and have no analogs in the anharmonic oscillator problem. Tadpole vertex Vtad 
(Jacobian source) is marked by (unhlled) open bullet. The signs of contributions and symmetry 
factors are indicated. 
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